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ABSTRACT 

We report the discovery of a correlation among three prompt emission properties of 
Gamma Ray Bursts. These are the isotropic peak luminosity L{ so , the peak energy 
(in vLv) of the time-integrated prompt emission spectrum -E" p k, and the "high signal" 
timescale To. 45, previously used to characterize the variability behavior of bursts. In the 
rest frame of the source the found correlation reads Li so oc #pk 2r G~45 49 - We find other 
strong correlations, but at the cost of increasing the number of variables, involving the 
variability and the isotropic energy of the prompt emission. With respect to the other 
tight correlations found in Gamma Ray Bursts (i.e. between the collimation corrected 
energy Ey and -E p k, the so-called Ghirlanda correlation, and the phenomenological 
correlation among the isotropic emitted energy Ei so , £" p k and the jet break time tbreak), 
the newly found correlation does not require any information from the afterglow phase 
of the bursts, nor any model-dependent assumption. In the popular scenario in which 
we are receiving beamed radiation originating in a fireball pointing at us, the found 
correlation preserves its form in the comoving frame. This helps to explain the small 
scatter of the correlation, and underlines the role of the local brightness (i.e. the 
brightness of the visible fraction of the fireball surface) . This correlation has been found 
with a relatively small number of objects, and it is hard to establish if any selection 
bias affects it. Its connection with the prompt local brightness is promising, but a 
solid physical understanding is still to be found. Despite all that, we find that some 
properties of the correlation, which we discuss, support its true existence, and this has 
important implications for the Gamma Ray Burst physics. Furthermore, it is possible 
to use such correlation as an accurate redshift estimator, and, more importantly, its 
tightness will allow us to use it as a tool to constrain the cosmological parameters. 
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1 INTRODUCTION 

A major breakthrough in the understanding of Gamma 
Ray Burst (GRB) has been the discovery of X-ray, optical 
and radio c ounte r parts in some long-dura t ion bursts 
JCosta et alJ Il997l: Ivan Paradiis et al. lTl997l iFrail et alJ 
I1997D . This allowed to measure their redshifts and real- 
ize that they are at cosmological distances. Since then, 
two potential and appealing possibilities were cherished: 
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2005). Both items face the problem of small number 
of events with measured z, since to determine it, deep 
optical/IR or X-ray spectra are needed. Regarding the 
latter item (ii), the large dispersion of the GR B energetics 
makes GRBs all but standard candles jFrail et al] boOll : 
iBloom. Frail fc Kulkarnl 12003ft . However, major advances 
in overcoming both difficulties can be done on the base of 
tight relations that connect GRB intrinsic energetics and/or 
luminosities with observed quantities. These relations serve 
as redshift indicators or, when z is known independently, to 
make GRBs standard candles for cosmographic purposes. 
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The first proposed GRB luminosity (redshift) indicators 
were based on the discovery of two empirical relations be- 
tween the burst isotropic-equivalent luminosity (Li so ) and 
the variability (V), which is a measure of the "spikiness" 
of the 7-ray light cu rve jFenimore fc Ramirez-Ruizl l200ot 
iReichart et ail 1200 J), and between Lj ao a nd the spectral 
lag. (r iai^ lNorris - Maxani. fc Bonnelll 12000'). Later. lAtteia I 
( 2003) constructed another redshift indicator based on 7-ray 
data alone linking the spectral parameters and the duration 
(T90) of the prompt emission, and he showed that it provides 
pseudo-redshifts accurate within a factor of two. In a study 
of the e nergetics of ORBs lElovd-Ronning fc Ramirez- Ruiz I 
(2002a) found the evidence that GRBs with highly vari- 
able light curves have greater uF v spectral peak energies. 
The existence of such a correlation and of the variability- 
luminosity correlation also implied that the rest frame GRB 
peak energy - E aea k is correlated w ith the intrinsic luminosity 
of the burst. lAmati et ail l)2002|) . by analyzing the spectra 
of BeppoSAX GRBs, found that the isotropic-equivalent 
energy radiated during the prompt phase (Ei so ) is corre- 
lated with the rest-frame peak energy of the 7-ray spectrum 
[E p k = £"pk S (l + z), E°^ s is in the observer frame]. However, 
the large scatter of this correlati on makes it difficult to use 
it as a reliable redshift indicator. lYonetoku et all l)2004l see 
also Ghirlanda et al. 2005) found that there is also some 
correlation between Li so and -Epk- 

Using the above mentioned empirical correlations, hun- 
dreds of GRB pseudo-redshifts were estimated from prompt 
7-ray observables and within a given cosmology. However, 
due to the lack of enough low-z (calibrator) GRBs, these 
correlations are themselves cosmology dependent (but see 
Ghirlanda et al. 2005a). This introduces a circularity prob- 
lem if the goal is to use these relations as distance-indicators 
and transform GRBs with measured z into cosmological 
rulers: the parameters of the relations depend on the cosmo- 
logical parameters that we p reten d to constrain. To avoid 
this circular logic, ISchaeferl i2003fl proposed to determine 
both the GRB and cosmological parameters by fitting si- 
multaneously the model relations to the data, and the cos- 
mological model to the constructed Hubble diagram (lumi- 
nosity distance dt, vs z). However, as he showed, an accurate 
determination of the cosmological parameters was not pos- 
sible due to the large dispersion around the Li so — V and 

- nag relations. 

The hope to use GRBs as cosmological rulers renewed 
after the finding by Ghirlanda, Ghisellini & Lazzati (2004a, 
hereafter GGL2004) that the Amati correlation becomes 
much tighter if one corrects Ei so for the collimation factor of 
the jet opening angle. The presence of an achromatic break 
in the GRB afterglow light curves is a strong evidence that 
the GRB emiss ion is collimated into a cone of semiapertur e 
angle 0j fe.g.. lRhoads1ll997l:ISari. Piran fc Halpern Ill99flft . 
where 6j is the rest frame time of the achromatic break in the 
lightcurve of the afterglow. Thus, the collimated-corrected 
7-ray energetic is i? 7 = (1 — cos 8j)Ei so . The _E 7 -_E p k corre- 
lation, determined with less than 20 GRBs, has been used 
to p ut constraints on some c osmological paramet ers avoid- 
ing IIGhirlanda et alJ l2004 bt iFirmani et alJ I2005T) and not 
avoiding foai. Liang fc Xull2004h t he circ ularity problem. 

Mor e recently iLiang fc Zhand l|2005l hereafter LZ2005) 
(see also lXu Il2005h investigated the correlation among iJiso, 
-Epk, and tbrcak without imposing any theoretical model and 



assumptions. They found a purely empirical tight correla- 
tion among these quantities which apparently suffers less 
uncertainties than the Ghirlanda relation, though both cor- 
relations are mutually consistent (Nava et al. 2006). These 
authors used the -E'iso--E p i c -tbreak correlation as a luminosity 
distance indicator and applied their own methods to avoid 
the circularity problem to constrain some dynamical and 
kinematical cosmolog ical parameters. 

As discussed in iGhisellini et alJ J20051) and LZ2005, 
the future of GRB cos mology is promising (but see 
iFriedman fc Bloom lEooih . the identification of more and 
more GRBs with accurately measured z, E°^ B and t^rcak 
being crucial. Unfortunately, the determination of tbreak re- 
quires expensive follow-up campaigns involving large tele- 
scopes. Up to now, fbreak has been measured for about one 
third (~ 20) of the GRBs with measured z (~ 60). This 
motivates the need of developing an astronomic al program 
entire ly devoted to measure z, _E p k and ibrcak llLamb et alJ 
2005). On the other hand, the discovery of tight correlations 
among the GRB energetics and the prompt 7-ray observ- 
ables alone (by-passing the need of measuring tbrcak) is of 
primary interest. This is the main aim of this paper. 

In §2, we present the GRB sample used in this work as 
well as the sample selection criteria. Using this sample, in 
§3 we revisit the correlation found by LZ2005 among Ei so , 
E p k and tbrcak- We will later compare this correlation with 
the one we find using observables which can be extracted 
by the prompt emission only. In §4 we present our method 
to characterize the variability of the prompt emission. The 
search for a multi variable correlation among energetics and 
the prompt 7-ray observables is presented and analyzed in 
§5, where we derive the main result of our work, namely a 
very tight correlation between the peak luminosity of the 
prompt emission Liso, the peak spectral energy E p ^ and the 
"high signal" timescale To. 45. In §6 we discuss how this newly 
found correlation can be used as a rather accurate redshift 
indicator. The theoretical implications of our findings are 
discussed in §7. Finally, in §8 we draw our conclusions. 



2 THE SAMPLE 

The study of the correlations among prompt emission ob- 
servables requires a careful selection of the GRB sample in 
order to work with a reliable and homogeneous set of infor- 
mation. The basic selection criteria we adopt are the knowl- 
edge of: 

(i) the redshift z; 

(ii) the peak flux P and the fluence F, better if defined 
both in the same energy range; 

(iii) the vF v peak energy, E°^ B , and a fitting to the spec- 
trum able to provide a reasonable determination of the bolo- 
metric correction. Note that E°^ B is derived in all cases using 
the time integrated spectrum. 

(iv) T0.45 and V. When possible, these parameters have 
been calculated with our own code using the published light 
curves (see below). 

The time is the prompt "h igh-signal" timescale, intro- 
duced bv lReichart et al. I j200ll hereafter R2001), defined in 
the rest energy range 50-300 keV (see §4.1). For compara- 
tive purposes, the knowledge of the observed tbrlak is also 



A tight correlation among prompt emission observables of long GRBs 3 



GRB 


Instrument 


z 


Ref. 


rpohs 

90 


Ref. 


'robs 
45 


Ref. 


DrCclK 


Ref. 


970228 


SAX/WFC 


0.695 


1 


80 


27 


2±1 


33 






970828 


RXTE/ASM 


0.957 


2 


146.6 


27 


14±2 


33 


2.2±0.4 


35 


971214 


SAX/WFC 


3.42 


3 


35 


27 


7.3±0.2 









980703 


RXTE/ASM 


0.966 


i 


102.4 


27 


17.2±0.3 





3.4±0.5 


35 


990123 


SAX/WFC 


1.6 


5 


100 


27 


17.2±0.1 





2.0±0.5 


36 


990506 


BAT/PCA 


1.307 


6 


220 


27 


14.05±0.06 









990510 


SAX/WFC 


1.619 


7 


75 


27 


5.00±0.05 





1.6±0.2 


37 


990705 


SAX/WFC 


0.843 


8 


42 


27 


5±2 


33 


1.0±0.2 


35 


990712 


SAX/WFC 


0.43 


7 


20 


27 




33 


1.6±0.2 


38 


991216 


BAT/PCA 


1.02 


9 


24.9 


27 


3.78±0.02 





1.2±0.4 


35 
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Uly/KO/NE 


4.5 


10 


110.1 


27 


4±1 


33 
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Uly/KO/NE 


1.058 


11 
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27 


5.2±0.5 


34 
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BAY /wFr 


2.14 


12 








o 


1.5±0.2 


39 
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HETE 


3.2 


13 


50 


28 


14.0±1.5 





3.0±0.4 


10 


020813 


HETE 


1.25 


14 


90 


27 


17±1 





0.43±0.06 


41 
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1.01 


15 


13 


28 


0.81±0.05 
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1.98 


16 
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28 


26±2 
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17 


316 


28 
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12 
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HETE 


0.1685 


18 


33 


28 


5.0±0.1 





0.5±0.1 


43 
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77 


28 
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20 
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21 


0.45±0.02 
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HETE 
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25 


21 


4.47±0.01 
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15 
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22 


8.8 


29 
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30 


32 
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Table 1. T° bs , T° b |, T° r b c s ak are expressed in seconds. References: (0) This paper; (1) Djorgovski et al. 1997; (2) Djorgovski et al. 2001; 
(3) Kulkarni et al. 1998; (4) Djorgovski et al. 1998; (5) Hjorth et al. 1999; (6) Bloom et al. 2001; (7) Vreeswijk et al. 2001; (8) Amati et 
al. 2000; (9) Vreeswijk et al. 1999; (10) Andersen et al. 2000; (11) Price et al. 2002; (12) Amati et al. 2004; (13) Hjorth et al. 2003; (14) 
Barth et al. 2003; (15) Vreeswijk et al. 2003; (16) Greiner et al. 2003; (17) Rol et al. 2003; (18) Greiner et al. 2003a; (19) Weidinger et al. 
2003; (20) Wiersema et al. 2004; (21) |http://space.mit.ed u.HETE/Bursts (22) Foley et al. 2005; (23) Tinney et al. 1998; (24) Galama 
et al. 1999; (25) Djorgovski et al. 2001a; (26) Bersier et al. 2003; (27) Ghirlanda et al. 2004a; (28) Sakamoto et al. 2005; (29) Blustin 
et al. 2005; (30) Jimenez et al. 2001; (31) Amati et al. 2002; (32) Mereghetti & Gotz 2003; (33) Reichart et al. 2001; (34) Guidorzi et 
al. 2005; (35) Bloom et al. 2003; (36) Kulkarni et al. 1999; (37) Israel et al. 1999; (38) Bjornsson et al. 2001; (39) Jakobsson et al. 2003; 
(40) Berger et al. 2002; (41) Klose et al. 2004; (42) Andersen et al. 2003; (43) Berger et al. 2003; (44) Jakobsson et al. 2004; (45) Stanek 
et al. 2005; (46) Bjornsson et al. 2001. 



desirable, although we did not require it as a selection crite- 
rion. For the few cases when the observed P is not available, 
we have estimated it from the light curve profile by using 
the published spectrum and fluence. 

Our sample includes mainly the GRBs reported in 
GGL2004 (15 GRBs); in some cases the data of these bursts 
have been updated with more recent published analysis. To 
the GGL2004 sample, we have added 7 GRBs. In Table 1 
we present the resulting sample of the 22 GRBs used in 
this paper, specifying the mission/instrument from which 
the data were acquired and the spectroscopically measured 
redshift z . For each burst we report the duration Tgo , the 
T0.4! timescale and the jet break time T^ak i n the observer 
frame. In most cases Tg 4I was computed from the GRB light 
curve as discussed below (§4.1). The data from R2001 where 
adapted to the 50-300 keV band by interpolation. We have 
compared our results with the ones of R2001 for the com- 
mon GRBs finding a good agreement. This comparison has 
been helpful for us to estimate the uncertainty of the data 
derived from R2001. 

In Table 2, we list the bolometric corrected fluence (F) 
and peak flux (P) and the corresponding energy range where 



they were computed. In Table 2 we also report the spectral 
parameter s a and for t hose b ursts whose spectrum was fit- 
ted by the lBand et al. I <ll993t) spectral model and only the 
a spectral parameter when the spectrum was fitted with a 
cutoff-powerlaw model. The fluence and the peak flux are 
given respectively either in units of erg cm -2 and erg cm -2 
s _1 , in round brackets, or in units of photon cm -2 and pho- 
tons cm -2 s _1 , in square brackets. This choice of homogene- 
ity, as well as the one of the same energy range for the fluence 
and the peak flux, have been made with the aim to reduce 
the effects of the spectral uncertainty on the Li so /Ei ao ratio. 

In Table 3 we report the rest frame -E p k, Ei so , Li so /Ei so 
ratio and V. The sample has been divided into two classes. 
Class 1 (16 events) comprises those GRBs with a good de- 
termination of the bolometric correction. Class 2 (5 events) 
comprises those events with uncertain bolometric correc- 
tions. In fact, the uncertain high-energy spectrum of Class 2 
bursts does not allow to constrain the high-energy Band 
spectral index (/3). In these cases spectral fits were per- 
formed either freezing (5 to a fiducial value (f3 = —2.3) or 
the fit was performed using a cutoff power law model. For 
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Figure 1. Isotropic emitted energy -Ej so as a function of 



,il.74 



-o.sr. 



peak break 



(Eq. 1.). The open circles represent the 15 bursts 
(i.e. all those with measured (break in Tab. 1 except 050525) used 
for the fit (solid line). The solid filled region represents the lcr 
uncertainty on the best fit (see Eq. 1). The fit is performed in the 
centroid defined by the data points (solid circle) where the er- 
rors on the best fit parameters are uncorrelated. The color codes 
are: blue for class 1 bursts, green for class 2 and brown for class 
3. For GRB 050525 we also report the values obtained assuming 
P = -2.1 and f3 = -5 (dotted lines). 



GRB 050525 we assign a special class, 3, because the deter- 
mination of (5 in this case is highly uncertain. 

Regarding the uncertainties, for those quantities taken 
from the observations (z, F, P, and E°]^ B ), we use the errors 
reported in the literature. For the derived quantities Ei so , 
Li so , and <& (see below), we calculate their uncertainties by 
propagating the errors under the assumption that they are 
uncorrelated. 

The conventional ACDM cosmology with Qm = 0.3, 
Qa — 0.7 and h = 0.7 is assumed here for the calculation of 
luminosity distances. 



3 THE £i so -£p k -£brcak CORRELATION 

As stated in the Introduction, our aim is to find a statis- 
tically significant (multi)variable correlation among GRB 
prompt emission observables. Previous to this, we present 
here for our sample the Ei so , E pk and t break correlation 
(LZ2005), and discuss some of its implications. We will also 
attempt later to compare and connect our results to this 
correlation. Using the same fitting procedure as in Nava et 
al. (2006) for those 15 events for which fbroak is known, we 



obtain the following best fit, in the GRB rest frame: 

1.74±0.10 



Sis, 



10 



53.11±0.04 



^break 

,10-°- 29 d 



( gj* v 

U0 248 keV/ 

-0.85±0.15 

erg, 



(1) 



with Xr=0-67 (for 12 dof). Notice that, differently from 
LZ2005, we do take into account for the fit the uncertain- 
ties in the data. The data and Eq. Q are shown in a 
two dimensional plot in Fig. Q In our sample, which has 
been constructed to find correlations among prompt vari- 
ables alone, GRB 020405 and GRB 021004 are not included 
because for them the high-energy spectral parameter /3 is 
larger th an -2 (GRB 020 405) or unknown (GRB 021004), 
(but see iNava etall l200rj) . We did not take into account 
GRB 050525 for the fit of Eq. due the large uncertainty 
in the parameter, but we plot it in Fig. the dotted lines 
correspond to values of f3 between —2.1 and —5.0. 

The multi variable correlation given in Eq. Q not only 
establishes a link between the 7~ray prompt and the af- 
terglow, but it may carry valuable information on the jet 
opening angle 8j. In the context of the fireball uniform jet 
model, this angle (supposed to be the same in the prompt 
and in the afterglow) is related to ibreak, the achromatic 
break time of the afterglow light curve (iRhoads I Il997t 
ISari. Piran fc Halpernl l999) . Depending on the circumstel- 
lar medium distribution, two models have been proposed. 
For the homogeneous medium (HM) with number density 



0j = 0.161 (t broak ) 3/8 



1/8 



(HM). 



(2) 



For a wind medium (WM) with a density profile n(r) = 
5 x 10 11 J 4r~ 2 g cm -2 (here A = 1 corresponds to the mass 
loss rate due to a wind of M w = 10 M© yr _1 and a wind 



velocity i\» 



lO' 5 km s" 



we have 



: 0.202 (tbreak) 1/4 



1/4 



(WM), 



(3) 



where 8j is expressed in radians and 7y 7 is the efficiency 
of the emission, i.e. the fraction of the kinetic energy of 
the fireball which is emitted in the prompt 7-ray phase. 
Such estimate of 8j implies the knowledge of z and of the 
luminosity distance. 

Interestingly enough, by introducing the found multi 
variable correlation Eq. Q in Eqs. (5) and (J^J, one may 
express 8j in terms of prompt observables. Thus, we obtain 
approximately that: 



171 3/2 

-bpk ; 

-Eiso 



(HM); 



a 2 -^pk 



(WM). 



(4) 



Defining the collimation corrected energy approximately as 
E y =E iso 8j, Eqs. © lead to 



E 1 oc £ pk 3/2 (HM); E-, <x E pk (WM). 



(•>) 



These relat i ons we re extensively discussed by GGL2004 and 
INava et ail teOOrt) . Finally, if we identify E^/E^ with the 
total number of photons N 1 emi tted by the prom pt, then 
an interesting conclusion follows ijNava et al.ll2006T) : 



1/2 



(HM); 



: const ~ 10 57 (WM). (6) 



In §7 we will use these results to compare the i?i S o-S p k-ibroak 
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[0.4+0.1] 


50-300 


17 
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-1.88±0.07 


-2.48i0.56 


(6.5±0.3)e-6 
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(1.8±0.1)c-5 


2-400 


5 


[40+4] 


2-400 


5 


031203 
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(2.4±0.2)c-7 
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18 



Table 2. a and j3 represent the photon spectral index of the model spectrum. The fiuence is given in units of erg cm" 2 (round brackets) 
or in units of photon cm -2 (square brackets). The peak flux is given in units of erg cm -2 s — 1 (round brackets) or photon cm -2 s _1 
(square brackets). The energy range of the fiuence and peak flux is expressed in units of keV. References: (0) This paper; (1) Ghirlanda 
et al. 2004a and refer ences therein; (2) Jimenez et al. 2001 ; (3) Amati et al. 2002; (4) Amati et al. 2004; (5) Sakamoto et al. 2005; (6) 
Atteia et al. 2005; (7) |http://space.mit.edu.HETE/Bursts| (8) Golenetskii et al. 2004; (9) Blustin et al. 2005; (10) Yonetoku et al. 2004; 
(11) The BATSE catalogue (Pacicsas et al. 1999); (12) Prontera et al. 2001; (13) Andersen et al. 2000; (14) Price et al. 2002; (15) Piro 
et al. 2005; (16) Yamazaki et al. 2003; (17) Paciesas et al. 1999; (18) Sazonov et al. 2004. 



or Ghirlanda correlation with the new one we will introduce 
in §5. 



4 CHARACTERIZING THE VARIABILITY OF 
THE PROMPT EMISSION 

The starting point, as mentioned in the Introduction, are 
the Liso-Spk or the Ei so -E p k correlations. Our aim is to see 
if the introduction of more prompt variables can produce 
significantly tighter correlations than these ones. Among the 
most relevant information from the 7-ray prompt emission 
is the lightcurve. 

The variability V, a quantity that estimates 
the "spikiness" of the ligh tcurve (see for details 
iFenimore fc Ramirez-Rurzl l2000t) . offers a way to par- 
tially quantify the information contained in the light 
curves. 

Here we introduce another variability indicator making 
use of the fiuence F to the peak flux P ratio. Such ratio gives 
a time, which compared with a burst timescale provides a 
further characterization of the lightcurve pattern variability. 
We calculate the F/P ratio using the bolometric corrected 
quantities (in the range of 1 — 10 4 keV). Using the prompt 



"high-signal" timescale, T0.45 ( see next subsection) we de- 
fine, in the GRB rest frame 

where To.45 = T °l|/(1 + z). $ is a scalar with a clear mean- 
ing in the observer frame, the GRB rest frame, and the fire- 
ball comoving frame. If one multiplies the numerator and 
denominator of the right term by 9?/E v k, then "I> reveals 
itself as the ratio of the number of photons emitted dur- 
ing the "high-signal" regime to the total number of emitted 
photons, which is particularly interesting in the comoving 
frame. A further advantage is that $ may be derived from 
observables easy to handle. 

4.1 Definition of variability 

The definition of the variability V requires the estimate of 
a smoothing time scale, T0.45, that, following R2001, we as- 
sume to be the time spanned by the brightest 45% of the 
total counts above the background. In practice, we estimate 
from the light curve the fraction of counts and the total time 
when the signal exceeds a given threshold. We identify the 
threshold when this fraction is 0.45, and we calculate the to- 
tal time during which the signal exceeds this threshold. The 
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Table 3. The peak energy E p ^ and the bolometric isotropic equivalent energy Ei ao are expressed in units of keV and erg, respectively. 
Ei BO and L; so are calculated in the 1-10 4 keV rest frame energy band. The Variability (Var) is calculated as described in the text (§4). 
References: (0) This paper (1) Ghirlanda et al 2004a (references therein); (2) Jimenez et al. 2001; (3) Amati et al. 2002; (4) Amati et 
al. 2004; (5) Atteia et al. 2005; (6) Barraud et al. 2003; (7) Sakamoto ct al. 2005; (8) Reichart et al. 2001; (9) Guidorzi et al. 2005; (10) 
Golenetskii et al. 2004; (11) |http://space.mit.edu.HETE/Bursts (12) Blustin ct al. 2005; (13) Sazonov et al. 2004; 



assumed energy range is 50 — 300 keV at rest. We use the 
recipe proposed by R2001 to pass from the observed energy 
range to the assumed rest one. Finally, we use the lightcurve 
time binning of HETE-II, 164-ms. Though for most of the 
BATSE GRBs the time bin is 64-ms, we preferred to take 
the larger HETE-II time bin for uniformity. In Table 1 we 
report the T0.4! calculated by us from the publicly avail- 
able lightcurves. In the few cases when the lightcurves were 
not available, we have used the time scales Tg\% reported 
in R2 001 and, for one event (see below), in lGuidorzi et ail 
l|2005h . 

The V parameter decreases systematically with the S /N 
ratio. In order to handle a quantity characterizing the signal 
it is necessary to disentangle the noise contribution. To this 
purpose R2001 suggested an analytic approximation. We in- 
troduce here an alternative method, with the aim to reduce 
the uncertainty on the variability noise correction avoiding 
any approximation. Our method is especially useful for data 
with a low S/N ratio. Based on the assumption that the 
photon shot noise obeys a Poissonian distribution, we use a 
Monte Carlo approach to simulate the noise contribution to 
the variability. 

The basics of our approach is as follows. Given a noise- 
less signal (zero order) with an intrinsic variability, we can 
introduce a first order realization of the signal by distorting 
it with its own Poisson noise. The first order variability will 
be different from the intrinsic variability. Then we construct 



a second order realization of the signal distorting the previ- 
ous first order realization with its own Poisson noise. This 
second order realization will have a new variability. Higher 
order realizations may be carried out up to a fixed number. 
The sequence of realizations may be randomly repeated and 
the average behavior and scatter of the variability sequences 
may be obtained. A simple statistical analysis allows to re- 
cover the intrinsic variability with some uncertainty from 
the higher (first, second,...) order realizations. Numerical ex- 
periments up to the third order realizations have been car- 
ried out successfully on artificial signals. We have found that 
the simple linear extrapolation Vo = 2Vi — V2 is enough to 
obtain a reasonable zero order estimate (here the index iden- 
tifies the order). In the case of a real signal, the zero order 
signal is unknown while its intrinsic variability is the result 
we are looking for. Actually we have just a first order real- 
ization with its intrinsic Poisson noise. Starting from here we 
apply the previous scheme building several times the second 
order realizations. The final result is the intrinsic variability 
with its uncertainty. 

We have taken some care to test our method using a 
real GRB light curve. To this aim, we have assumed that the 
observed light curve is a pure zero order signal (i.e. ignor any 
Poisson noise and call Vo the variability of this light curve). 
Then we have distorted this signal once with a Poisson noise 
obtaining a first order realisation (Vi), and finally on this 
one we have applied our method (i.e. derive V2 by a Monte 
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Carlo method, extrapolate back and estimate the average 
Vo and its error). The result may be compared now with the 
variability of the "zero" order signal Vo- In the case of the 
rather noisy GRB 971214 light curve, Vo = 0.1169 lies inside 
the obtained V = 0.116±0.004 (la error), while for the high 
signal GRB 991216 light curve we have V = 0.1438, which 
is consistent with Vo = 0.1445 ± 0.0010. 

In principle our variability estimate follows a different 
procedure from the original method proposed by R2001, 
therefore a comparison between the two methods is desir- 
able. We have used the GRBs of our sample in common 
with the R2001 sample to compare the values of T0.45 and 
V. From the results given in Table 1 and Table 3 it can 
be seen that our method gives T0.4I and V in general good 
agreement with those published by R2001. The only differ- 
ence is a smaller uncertainty in our estimate of V. The fact 
that our results agree with R2001 justifies the inclusion in 
our sample of four GRBs whose V and T0.45 have been esti- 
mated by R2001 only 

In a recent paper. IGuidorzi et alJ i2005l) have presented 
new calculations of V for a sample of 32 GRBs. For some 
GRBs in common with our sample, the processed lightcurves 
used by these authors could differ from the ones used by 
us. Besides, for the BeppoSAX/WFC events they apply two 
further "instrumental" corrections in the calculation of V 
that we did not apply because of the lack of information. 
According to the authors, these corrections affect V signif- 
icantly only for relatively short GRBs exhibiting sharp in- 
tense pulses. On the other hand, our calculation of V and 
its uncertainty includes a new method to treat the noise 
correction that improv es the results especiall y for data with 
low S/N ratios, while IGuidorzi et alJ (120051) use the same 
analytical approach of R2001. We h ave compared Tfi\% an d 
V for the 5 GRBs in common with IGuidorzi et al.l {2005). 
Some discrepancies have been found. Therefore, for reasons 
of homogeneity, we hav e included in our sam ple the mea- 
sured T0.4! and V from IGuidorzi et all J2005li only for one 
event, GRB 000911. 



5 A NEW TIGHT CORRELATION 

Our aim is to find correlations among prompt 7— ray quanti- 
ties alone, with the smallest scatter. For example, we would 
like to explore the possibility to improve the Li so -E p k (Ei so - 
E p k) "power" ("energy") correlations by introducing the 
prompt emission variables To. 45, V, and <fr, which partially 
characterize the 7-ray light curve (see §4). We remark that 
in our statistical analysis we take into account the uncer- 
tainties of the data. Notice that the variables always refer 
to the GRB rest frame. Our main criterion to judge about 
the goodness of the fit will be the value of the reduced x 2 ■ 
Among other combinations, we look for a multi vari- 
able linear correlation of Log$ as a function of a combi- 
nation of LogiSpk, Logi5i so , LogXb.45 LogV (Liso is already 
contained in the variable $). This approach also allows us to 
explore whether the lightcurve variables, $ and V, are cor- 
related among them or not. Only the 15 high-quality (class 
1) events with accurately measured variability V are used 
in the regression analysis. The best fit is: 

1.51±0.25 




Log(EJjf' E£"> T°;« V°-«) 



Figure 2. Multi variable linear regression analysis based on the 
15 GRBs of class 1 (open circles). The best fit (solid line) corre- 
sponds to Eq. |S] and the shaded region is the la uncertainty of 
this equation. The latter is computed in the barycenter of the data 
points (solid black point). Color codes have the same meaning as 
in Fig. El 
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E, 



pk 



10 2.38 keV 



with Xr=0.62 (for 10 dof). Fig. H shows the data and the 
best fit given by Eq. @ in a two dimensional diagram. The 
uncertainties in the data, showed with error bars, were taken 
into account in the regression analysis. The (yellow) shaded 
area corresponds to the la confidence interval of the regres- 
sion line (Eq|HJ. The data used for the fitting are marked 
with open circles. The class 2 events and GRB 050525 are 
also plotted in Fig. |21 although they were not used in the 
fit. 

From Eq. © one sees that $ is correlated with V 
very weakly, since its exponent is almost consistent with 
zero. We also note that _Bi so is not an important variable in 
the correlation represented by Eq. (JHJ, because it appears 
with almost the same power on both sides of this equation. 
Therefore, the significant variables are Li so , E p k and Tocs- 
in a first approximation, from Eq. (SJ we may infer that 
Li so oc _E pk ■ Xb.45 - ' • 

We therefore proceed to perform the multi variable re- 
gression analysis (taking into account the uncertainties) for 
LogLiso as a function of LogL p k and LogTb.45. For the same 
events used previously, we obtain a fit with Xr=0.67 (for 
13 dof) that, compared with the Xr— 0.62 of the correlation 
given by Eq. JHJ, implies again that the variables V and _Ei so 
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Figure 3. Multi variable linear regression analysis based on 19 GRBs (class 1 [blue crosses] and class 2 [green crosses] events in Tab. 
3). The best fit, corresponding to Eq. [5] and its lcr uncertainty (computed in the barycenter of the data points - solid black point) are 
represented by the solid line and by the shaded (yellow) region, respectively, using the same convention of Fig. HI 



are redundant for the correlation in Eq. |H| In order to use 
as many data as possible, we also performed the regression 
analysis including 4 events of class 2 (GRB 030429 lacks a 
measure of Tq 4I), for a total of 19 events. The best fit is: 

/-rp s 1.62±0.08 

r _ 1n 52.11±0.03 / -frpk \ 

LiS ° ~ 10 W-37 keV J 
/ rp . -0.49±0.07 



with Xr=0.70 (for 16 dof). Fig.EJshows the data, the best fit 
and its la confidence interval (shaded area) in a two dimen- 
sional diagram. GRB 050525 was not used for the regression 
analysis but it is plotted in Fig.[!|] 

The correlation Z/i S o--E p k-To.45 is very tight, which can 
be appreciated from Fig.0] where we show the scatter (taken 
as the orthogonal distance of the data points to the best 
fit line). Fitting the scatter distribution with a Gaussian 
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Figure 4. Histogram of the (logarithmic) scatter of the data 
points around the best fit of Ea.l9l flight gray) and Eq.Q (darker 
gray). We superpose the two Gaussian fits. 

yields a (logarithmic) dispersion a — 0.06, smaller than the 
corresponding dispersion of points around the correlation 
described by Eq.0(a Gaussian yields a — 0.08). 

The correlation of Eq. is a significant improvement 
(i.e. it h as a smaller sca t ter) o f the Li BO -E p k correlation 
found bv lYonetoku et all (I2004T) . In fact, for the updated 
version of the Yonetoku correlation, Ghirlanda et al. (2005) 
found a la dispersion of 0.25. We have explored other multi 
variable correlations and the one given by Eq. © was indeed 
the tightest with significant variables. 

The fact that Xr is smaller than 1 might be a conse- 
quence of an overestimate in the observational uncertainties 
of the data as well as of unaccounted correlation among 
the variable errors. Indeed, we have assumed that the errors 
in the variables are uncorrelated (§2). If the errors instead 
have some degree of correlation among them (which is very 
likely), then we are overestimating their contribution to Xr- 

Fig. El shows the main result of this paper, namely the 
existence of a very tight empirical relation among the peak 
luminosity of the GRB (Li so ) and the prompt 7-ray quanti- 
ties _E p k and To. 45. This relation is as tight as the Rso-Epk- 
tbreak one (see §3). 

In Fig. E|we display the available GRB light -curves fol- 
lowing the same sequence of data points (from top-right to 
bottom-left) plotted in Fig. [3] The different patterns appear 
mixed, and there is no hint of any sequence. 

We have also performed a linear regression analysis to 
find the best fit in our sample for Li so -£ , p k-T9o , i.e. using 
the variable T90 at rest instead of To. 45. The obtained mul- 
tivariable correlation is similar to the one given by Eq. ||UJ 
but with a significant larger scatter. 

5.1 The lowest— luminosity GRBs in the 

£iso-Epk-To.45 diagram 

There are four GRBs that were not included in our sam- 
ple: GRB 980425 (a rather peculiar burst), GRB 990712 
and GRB 010921 (they were not included in our sample 




Figure 6. Similar to Fig. |3 but including the "outliers" 
GRB 980425, GRB 990712 GRB 010921 and G RB 031203. If 
the X-ray E p ^ reported by I Watson et alj 120061) is used, then 
the position of GRB 031203 (lower limit calculated for E p k = 
200 ± 4QkeV) shifts close to the correlation given by Eq. |5] 



because they have a too uncertain determination of V"), 
and GRB 031203 (another peculiar burst, as GRB 980425). 
All these bursts are low luminosity events. It is interest- 
ing to see a posteriori where they lie in the Logii so vs 
Log^ 162 ^- ' 49 ) plot. 

Fig. |S| shows the same correlation of Fig. [3] but we 
include here these four low luminosity events. Curiously 
enough, they break the main correlation at its lowest end, 
below Z/i so ~ 10 51 ' 2 erg s -1 , giving rise to a kind of low lu- 
minosity branch (determined mainly by a minimum value of 
-Epk and, at a less ext ent, by a maximum v alue of T0.45). 

In a recent paper. Iwatson et alJ i2006l) have shown that 
GRB 031203 might have a soft X-ray spectral component, 
at ~ 2 keV. This component carries a fluence which is even 
higher than the 7-ray one. Therefore it is this soft compo- 
nent which can determine the peak energy -E p k. 

If we use i? a k of the X-ray component reported in 
IWatson et ail |2006), then the position of GRB 031203 in 
Fig. [(J shifts to the left as indicated by the arrow in Fig. |S| 
In this case GRB 031203 lies close to the multi variable cor- 
relation given by Eq. ©. Furtherm ore, the reported sp ectra 
of GRB 990712 and G RB 010921 llAmati et al] (120021) and 
ISakamoto et al. I (l2005ll . respectively), seems peculiar: in the 
case of GRB 990712 the spectrum is rather flat in vF(y), 
making the determination of SpkUncertain, and in the case 
of GRB 010921 the spectrum has some local bumps at en- 
ergies lower than the reported £" p k (Barraud et al. 2003). It 
is then possible that also these two sources lie very close to 
the main correlation shown in Fig. |S| 



10 Firmani et al. 




Time (rest frame) 



Figure 5. The sequence of the available GRB light curves of our sample. The order follows the distribution of GRBs along the sequence 
of Fig.EJfrom top-right to bottom— left. 



6 THE L iso -E pk -T 0A5 CORRELATION AS A 
REDSHIFT INDICATOR 

The tightness of the multi variable Liso-Epk-To^s correla- 
tion (see Fig. |1J encourages us to use it as new distance 
(redshift) indicator for GRBs with no measured redshifts. 
In fact, the great advantage of this correlation is that in- 
volves only prompt emission quantities which can be easily 
derived from the 7-ray lightcurve and spectra of GRBs with 
unknown z. 

This way, we could have a statistically sig- 
nificant GRB luminosity-redshift diagram to infer 
the GRB luminosity fu nction and the GRB f orma- 
ti on rate history (e.g.. ISchaef er. Dene fc Ban d! l200li 
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Llovd-Ronning. Fryer fc Ramirez^ui3]2002 
20041: lYonetoku et al.ll2004l) . 

The bolometric corrected peak flux, in the observer 
frame, is connected to the bolometric corrected isotropic lu- 
minosity Liso by: 



L[ so 



(10) 



where Dl(z) is the luminosity distance. The cosmologi- 



cal redshift and time dilation effect impose that _E p k = 
E°l s (l + z) and T0.45 = T °«/(l + z). However, the use 
of T0.45 needs a particular care. In fact, as mentioned in 
§2.1, To.45 (and consequently To. 45) is calculated from a light 
curve defined over an energy range which is fixed in the rest 
frame of the GRB (see also R2001). Using the approximation 
given by R2001 that takes into account the narrowing of the 
light curve temporal substructures at higher energies, we ob- 
tain TSab = 20.45/(1 + zf A , where f °^l is measured in the 
same energy range but defined now in the observer frame. 
It follows that Tb. 45 measured in the source rest frame is 
To.45 = T ( f 4 !/(l + z) - 6 . 

Once we know how the rest frame quantities Li so , E p k 
and To. 4 5 transform with the redshift z, we may express 
the I/iso-Spk-To.45 relation (Eq.|5J in terms of the observer 
frame equivalent quantities: 



P4ttDI(z) 



10" 



T° k bs (l + ^) 



T~obs 
0.45 



10 2 - 37 keV 

-0.49±0.07 



1.62±0.08 



(l + z)°- 6 10 -' 



erg s (11) 
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986 
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Table 4. Spcctroscopically measured redshifts (z) and pseudo— 
redshifts (z) derived according to the procedure described in 
Sec. 7 for 20 GRBs of our sample. Az/z represents the uncer- 
tainty of the derived pseudo-redshifts. 




Figure 7. The line shows the pseudoredshift given by Eq. 1151 
while the data show the observed redshifts for the 19 GRBs of 
the sample of Fig. 151 



We indicate the ratio between observer frame measured 
quantities as $ and we isolate the redshift dependences in 
the term /(«): 



Fobs 1 



(13) 



so that Eg. II II reduces to: 

f(z) = C9, (14) 

where C is fixed to 10 6 48 such that D L , P, E°£ s , and T ol f s 
are expressed in unities of Gpc, erg cm -2 s _1 , keV and sec, 
respectively. 

For the ACDM cosmological model assumed in this pa- 
per, the inversion of Eq. 1 IH allows us to estimate the 
pseudo-redshift z of a given GRB whose P, E°[f , and To.« 
have been measured from the prompt emission 7-ray light 
curve solely: 



-z = r 1 {c^). 



(15) 



/(*) = 



(1- 



(12) 



In Fig.Qwe plot Eq. 1141 as Log^ vs. Log(z) (solid line). 
The sample of GRBs used to derive the Liso-.Epk-T0.45 cor- 
relation are also plotted (stars). The notably small scatter 
of the data points around the curve is a proof of the reliabil- 
ity of the redshift indicator inferred from the Ei so -Epk-Tb.45 
correlation. In Table 5 we give for the 19 GRBs of our sam- 
ple the spectroscopically measured z, the pseudo-redshifts 
z derived from Eq. 1151 and the relative uncertainty between 
these two quantities. We note that the 68% of our pseudo 
redshifts (i.e. 13 out of 19 sources) are within the 25% of the 
real value. In Fig.|H]we compare the distribution of the Az/z 
values obtained with Eq. 1151 with the one found by Atteia 
(2003). One can see that the our distribution is somewhat 
narrower, and especially that the central value is closer to 
zero. 

Moreover, as can be inferred from the data reported 
in Table 4, the relative uncertainty on the derived pseudo- 
redshifts does not seem to be correlated with the redshifts. 
This suggests that the Liso-Epk-Tb.45 correlation can be 
used to estimate either low and high values of z. 

Finally we note that the Tq 4I parameter can be derived 
from the prompt emission lightcurve straightforwardly. On 
the other hand, the use of Eq. 1151 to estimate z requires the 
measure of the GRB peak spectral energy E°^ B from the fit 
of a broad-band spectrum. 



7 THEORETICAL IMPLICATIONS AND 
CHALLENGES 

7.1 Comparison with other correlations 

The most important correlation found in this paper among 
Li so , .Epk and To. 45 can be compared with the other tight 
correlations between the emitted energy or luminosity, i? p k 
and ibreak- Consider first the Amati relation between i?i so 
and Epk- It has a relatively large scatter, which "collapses" 
when calculating the collimation corrected emitted energy 
E 1 . In this case, the "collapsing parameter" going from the 
Amati to the Ghirlanda relation is the jet opening angle, 
measured by ibreak- If indeed 0j is well measured by tbreak, 
then the existence of the Ghirlanda relation implies the exis- 
tence of the phenomenological i?iso-E p k-ibreak tight relation 
(Eq.0 (and vice-versa). 

Consider now the Yonetoku relation between Li so and 
-Epk. As detailed in Ghirlanda et al. (2005), it has a scatter 
which is comparable to the scatter of the Amati relation. 
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This paper 
Atteia (8003) 



-1 
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Figure 8. Histograms of the values of Az/z obtained with the 
method discussed in this paper (hatched, 19 GRBs) and by by 
Atteia (2003, dashed line, 17 GRBs). 



But if we now calculate L 7 (correcting for the jet opening 
angle), then we do not find a correlation as tight as the 
Ghirlanda one. Therefore the "collapsing parameter" which 
worked for the total energetics does not work for the peak 
luminosity. The correlation found in this paper demonstrates 
that a better "collapsing parameter" for the peak luminosity 
is Tb.45- 

From this we can already conclude that To. 45 is not 
equivalent to the jet opening angle: if it were, then both 
To. 45 and 8j (or ibreak) should work as collapsing parameters 
without "caring" if dealing with total prompt energetics or 
peak luminosities. 

Therefore the two tight empirical multi variable corre- 
lations Liso-Epk-To.45 (Eq. El and £'i so --Epk-ibrcak (Eq. EJ 
represent two different aspects of the GRB physics. 



of the -Eiao--E p k-tbreak an d the Z/iso-Spk-To.45 relations. We 
can divide both sides of Eq. |§]by $ (see Eq. [SJ so that E[ BO 
becomes the dependent variable in Eq. |5] as it is in Eq. 
A more rigorous result can be obtained by fitting (with the 
multi variable regression technique) Log£i so as a function 
of Log_E p k LogTo.45 and Log<3>. The best fit, obtained for 
our sample (including class 1 and 2 GRBs for a total of 19 
events), is: 



E iso = 10 



52.90±0.03 



( g P k y 

U0 2 - 36 keV/ 

0.51±0.07 



10°«sJ \ 



io- 



-1.00±0.15 



erg (16) 



with Xr= 0.64 (for 15 dof). This multi variable correlation is 
completely equivalent to that of Eq. @ . By eliminating _E p k 
in Eqs. I|16|l and we can express ibreak as a function of 
of $ and To. 45 (the dependence on i?i so becomes negligible): 

( $ 2 \°' 6 

ibreak ~ 4 — days (17) 

\-t0.45 / 

Ea . 1 1 71 establishes a clear connection between prompt emis- 
sion local quantities (<£> and T0.45), which are related to the 
7-ray lightcurve variability, and an afterglow emission global 
quantity, i.e. ibreak- 

LZ2005 suggested that ibreak could be related to the 
prompt (7-ray) emission rather than to the afterglow one. 
ibreak enters in the Ghirlanda correlation by means of the jet 
opening angle (see §3), and, as discussed in LZ2005, it raises 
the question of how a global quantity (i.e. 8j) conspires with 
Eiso to affect -E p k, which, instead, are both local quantities. 

It could be that ibreak is determined by only prompt 
emission local quantities. This i s as a lso suggested by the 
reported iSalmonson fc Galamal (|2002)) correlation between 
the pulse lag of the prompt 7-ray emission and ibreak. This 
might justify the existence of an empirical correlation be- 
tween the jet angle (derived from ibreak) and the local prop- 
erties Ei BO and _E p k. Therefore, our finding that ibreak is 
closely related to the prompt lightcurve quantities $ and 
To. 45 (Eq. 1171 could be the basis of the Ghirlanda and the 
i?i so -i5pk-tbreak correlations. 



7.2 "Local" versus "global" properties 

It is of some interest to make a distinction between local 
and global properties of GRBs. If the fireball has a semi- 
aperture angle Qj > 1/T, we can see only a portion of the 
emitting surface (or volume) . The peak flux, the fluence and 
the timescale To. 45 are defined as local quantities since they 
correspond only to the "observable" surface (or volume). 
This implies that also Li BO and Ei BO are basically local quan- 
tities. We can then introduce the concept of "local bright- 
ness" characterizing the flux emitted by the visible fraction 
of the fireball. On the other hand, the jet opening angle 
6j (and consequently ibreak) is a global quantity. This im- 
plies that i? 7 and the Ghirlanda relation are expressions of 
a global property of GRBs. On the other hand, the Lj so - 
i? p k-To.45 relation involves the local brightness (entering in 
the definition of Li so ), which is a local property. Further- 
more, Liso is an instantaneous quantity, not time integrated 
(contrary to Ei so or E^). 

We now discuss some implications of the combination 



7.3 The Li so — £pk— T0.45 correlation in the comoving 
frame 

Another important aspect concerns the transformation of 
the relation Liso-Epk-To.45 from the GRB rest frame to the 
fireball comoving frame. Due to blueshift, E pk and the co- 
moving (primed) i? pk are related by _E pk = E' pk 8, where 
S — l/[r(l — /3cos£)] is the relativistic Doppler factor, 
and £ is the angle between the line of sight and the ve- 
locity vector. Similarly, we have To. 45 = Tq A5 /S. For the 
isotropic equivalent luminosity we have two possible trans- 
formations, according if the fireball is collimated in an angle 
6j > 1/T ("standard" fireball) or if its relevant aperture 
angle 8j < 1/T [e.g. it is made by "bullets" (Heinz & Begel- 
man 1999; sub-jets (Toma, Yamazaki & Nakamura, 2005); 
or cannonballs (e.g. Dar & De Rujula 2004)]. In the former 
case we have Li so = <5 2 i( so , while in the latter Li so = S 4 L' iso . 
This is due to the fact that, for 0j > 1/r, radiation is colli- 
mated in a cone of semiaperture angle 6j, while, if 8j < 1/T, 
the collimation angle becomes 1/r. 
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Eiso 



For ease of discussion, let us approximate Eq. [5] as 
It is clear that, for 9j > 1/r, the 



rp3/2rp 

-^pk 1 0.45 



1/2 



same relation holds in the comoving frame: the 5 fac- 
tors cancel out. Instead, for 9j < 1/r, we have L[ so oc 



s.3/2 



(To, 



-1/2. 



there is a rather strong dependence 



<T 2 (E pk , 
on 8. 

In other words, Eq. [3] is "Lorentz invariant" only for 
"normal" fireballs (i.e. those with 9j > 1/r). In the oppo- 
site case (8j < 1/r) if the V factor varies more than a given 
(small) factor from event to event, then we have a conflict 
with the tightness of our relation. We can therefore conclude 
that the standard fireball scenario is favored with respect to 
models in which we see the entire emitting surface (cannon- 
balls, sub-jets and bullets, of angular extension 9 < 1/r). 

If the "Lorentz invariance" can be used as a discriminat- 
ing guide among the two possible forms of the Ghirlanda re- 
lation (favoring its "wind-like" version, which linearly links 
E 7 and E p k), then we have, from that relation, that the 
number of relevant photons of the prompt emission (those 
contributing the most to E 7 ) must be the same in dif- 
ferent bursts (see Eq. 6 and Nava et al. 2006). From the 
Liso-iSpk-ro.45 relation found in this paper, instead, we de- 
rive a relation between the local (and instantaneous) maxi- 
mum brightness of the fireball, the To. 45 timescale and E p k. 
This relation, in the comoving frame, reads 



(£pk) 3/2 



r n~ ~ r 



F' 



P k 



^0.45 



1/2 



(18) 



where r is the radius of the fireball when it emits the peak 
brightness F' and the photon peak flux n 7 . From this per- 
spective, we can see that the Ghirlanda relation and the 
relation found here (Eq. 9) are complementary, describing 
two different aspects (albeit hopefully related) of the GRB 
physics. 



8 CONCLUSIONS 

GRBs are extremely interesting objects: (i) their nature and 
physics are still a mystery to solve, and, (ii) GRBs, being 
the most powerful explosions in the Universe, can serve as 
a valuable cosmological tool. 

An important step forward in the study of GRBs has 
been recently achieved through the discovery of several cor- 
relations among their intrinsic properties. 

Through the analysis of a sample of 22 GRBs, we have 
discovered a new very tight correlation among 7-ray prompt 
emission properties alone (Sec. 5). The rest frame prompt- 
emission quantities that we have correlated are: the bolo- 
metric corrected energy and luminosity Ei so and Lj ao , the 
prompt emission spectrum peak energy -Epk, the variability 
of the 7-ray lightcurve V , and the "high signal" timescale 
To. 45 (as defined in Sec. 4.1). We have adopted in this pa- 
per an improved method to estimate V and To. 45 from the 
prompt emission lightcurve. Moreover, we adopted an inte- 
gral measure of the light curve variability as defined through 
the ratio <!> = Ej so Tb.45/Ei S o- The number of objects avail- 
able to find these correlation is unfortunately still small, and 
a robust physical explanation is still to be found. Bearing 
that in mind, the main conclusions of our study are: 

• Long GRBs with measured redshifts obey a very tight 



rest frame multi variable correlation which involves three 
prompt emission quantities: Ei so oc -Ep k 62 T - 4 5 49 (Eq. [5] and 
Fig. 3). The scatter of this correlation is comparable to that 
of the -Eiso, £ P k and tbrcak (LZ2005; Fig. 1). Nonetheless, 
the newly found Ei so -E p k-To.45 correlation involves only 
prompt-emission quantities. 

• The four GRBs (980425, 990712, 010921 and 031203) 
which populate the low-luminosity tail of the Ei so distri- 
bution of the bursts of our sample, were not included in 
the derivation of the above correlation. This is because of 
their arguable nature and/or high uncertainties in their 
observational parameters. However, with respect to the 
Li so -Epk-ro.45 correlation, they define a low luminosity 
"branch" (Fig. 6). Whether the origin of this branch is phys- 
ical or merely a consequence of the lack of spectral data at 
low energies (as for GRB 031203), is an intriguing question 
which requires further analysis. 

• Other tight correlations among a larger number of 
prompt emission quantities were also found by our multi 
variable analysis. For example, we reported a correlation of 
$ as a function of E pk , E iso , T 0A5 and V (Xr= 0.62; Eq.|U 
and Fig. 2), and of Ei so as a function of E p k, To. 45 and $ 
(Xr= 0.63; Eq. 1161 . However, these extra variables (with 
respect to the Ti so -E p k-To.45 correlation), only slightly im- 
prove the best fit Xr, at the cost to include a larger number 
of variables. 

• The tightness of the multi variable correlation 
Liso-Epit-To.45 allows us to use it as a cosmology-dependent 
redshift indicator for GRBs. For the concordance cosmolog- 
ical model used here, we presented the formulas necessary 
to infer pseudo-redshifts given three observables: the bolo- 
metric corrected P, E°^ B , and T0.45. This method represents 
an improvement with respect to previous methods. 

• We have re-derived the Eiso-Ep^-tbrcak relation 
(LZ2005) for the 15 GRBs out of our sample with known 
tbrcak- Then, we connected this relation to one equivalent 
to the Tiso-Epk-To.45 relationship (E iso [E pk ,t brca-k ,&}) and 
found tbrcak ~ 4($ 2 /To. 45 ) 0,6 days. This dependence es- 
tablishes a connection between prompt local quantities ($ 
and To. 45, related to the 7-ray lightcurve variability), and a 
global quantity associated to the afterglow (tbrcak), suggest- 
ing that tbrcak (and 9j) could be determined mainly by the 
prompt emission local properties. 

• In the context of the standard fireball scenario, the 
Lis -Epk-To.45 correlation is roughly "Lorentz invariant". 
This allows us to go deeper and relate it to the local comov- 
ing fireball surface brightness associated to the fireball area 
that the observer sees. Therefore, while the Ghirlanda rela- 
tions seems to be related to a global feature of the fireball 
(the number of emitted photons), the Tiso-Epk-To.45 rela- 
tionship could be describing a local feature of the fireball 
related to its brightness. 

The results found in this paper reveal the existence of 
tight correlations among quantities related to the prompt 
emission of long GRBs. On one hand, these correlations will 
contribute to understand better the physics of GRBs. On the 
other hand, they can be used to infer pseudo-redshifts with 
good accuracy for hundreds of GRBs or, for those events 
with known 2, serve as "standard candles" for cosmographic 
purposes. The main challenge is the acquirement of obser- 
vational data related to the 7-ray prompt light curves and 
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spectral distributions. In particular, detector flux sensitivi- 
ties extending to energies as high as 500-1000 keV are cru- 
cial, also for revealing the effects of possible observational 
bias. 
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